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ABSTRACT 

We make predictions for the pairwise velocity dispersion (PVD) of galaxies with models that are constrained 
to match the projected correlation function and luminosity function of galaxies in the Two-Degree Field Galaxy 
Redshift Survey (2dFGRS). We use these data to constrain the halo occupation distribution (HOD) of 2dFGRS 
galaxies, then calculate the PVD by populating the halos of a large-volume, high-resolution N-body simulation. 
We examine the luminosity and scale dependence of the predicted PVD. At large scales the PVD is flat, but at 
r ~ 2 h Mpc the dispersion rapidly rises as the dominant source of galaxy pairs transitions from two-halo pairs 
to pairs within a single halo. At small and large scales, r < 1 /r'Mpc and r > 3 /r'Mpc, we find that the PVD 
decreases with increasing galaxy luminosity. This result is mostly driven by the fraction of satellite galaxies 
/sat. which is well-constrained by the shape and amplitude of the correlation function. We find / sat ~ 25% for 
galaxies fainter than L*, while for brighter galaxies the satellite fraction rapidly declines, creating the decrease 
in the PVD with luminosity. At r = 1 A -1 Mpc, the PVD has no dependence on luminosity because satellite 
galaxies dominate the statistics for both bright and faint objects. Recent measurements of the PVD in Fourier 
space using the "dispersion model" have reported a strong decline in PVD with increasing luminosity at k = 1 
Mpc -1 h. We test this method with our HOD models, finding that the dispersion model roughly recovers the 
shape of PVD with scale, but that there is no consistent comparison between the PVD at a given A:-scale and 
the true dispersion at a given value of r. This results in a luminosity dependence in £-space that is stronger than 
in configuration space. The luminosity dependence of the HOD results in Fourier space are consistent with 
those measured at k = 1 Mpc" 1 h; thus the recent measurements of the PVD are fully explainable in the context 
of halo occupation models. The normalization of the PVD is lower than predicted by our fiducial model, and 
reproducing it requires a lower value of Q m (~ 0.2 instead of 0.3), a lower value of as (~ 0.7 instead of 0.9), 
or a significant bias between the velocity dispersion of galaxies and dark matter in massive halos. 
Subject headings: cosmology: theory — galaxies:halos — large scale structure of the universe 



1. INTRODUCTION 

Galaxy peculiar velocities, their motions with respect to the 
smooth Hubble flow, offer a unique window into both cos- 
mology and galaxy formation. The galaxy velocity field re- 
flects the matter density of the universe, fL,, and the ampli- 
tude of matter fluct uations ( Peeble*sl ["l976l ISargent & Turner! 
119771 lKaiseJ[T98 7). This is true at large scales, where the 
infall velocities of galaxies toward overdense regions are de- 
termined by the overall amount of matter and its distribution 
at large scales, and it is also true at small scales, where the 
random motions of galaxies of galaxies are determined by the 
shape and normalization of the mass function of dark matter 
halos ( Sheth 1996). But peculiar velocities also provide infor- 
mation about galaxy formation and evolution; just as galaxies 
are a biased tracer of the spatial distribution of dark matter, 
they can also be biase d with respect to the velocity distri- 
bution of dark matter dBerlind et ai] 120031 lYoshikawa et alJ 
120031 Faltenba cher et alJl2005l) . Pairwise velocity statistics 
generally differ even if galaxies have the same peculiar veloci- 
ties as local dark matter particles because spatial bias changes 
the weighting of those velocities. This is of particular im- 
portance at small scales, reflecting the accretion histories and 
evolution of galaxies within groups and clusters. A full model 
of galaxy formation must reproduce both the real-space clus- 
tering of the observed galaxy distribution and the distribution 
of velocities, and the dependence of these distributions on lu- 

1 Kavli Institute for Cosmological Physics, University of Chicago 

2 Institute for Astronomy, University of Edinburgh 

3 Department of Astronomy, Ohio State University 

4 Theoretical Astrophysics, Los Alamos National Labs 



minosity and galaxy type. 

In this paper we make predictions for the pairwise veloc- 
ity dispersion of galaxies (hereafter PVD, or a v ) based on 
the observed real-space clustering of galaxies in the Two- 
Degre e Field Galaxy Redshift Survey (2dFGRS: IColless et alJ 
l20f)lT) . We use the Halo Occupation Distribution (HOD) 
modeling method, in which the relation between galaxies 
and m ass is described at the level of virialized dark matter 
halos (Kauffmann et al. 1997: |Jing et al-||1998: Benso net alJ 
| 2000t ISelia1dl200(l iPeacock & Smithll2000t iMa & Frvll20C 
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Scoccimarro et alJ 1200 It iBerlind & Weinberg! 120021) . The 
mass function and spatial clustering of dark matter ha- 
los themselves, for given cosmological parameters, are 
well-determine d from numerous numerical and analytic 
studies (e.g.. ISheth & Tormed 119991: jJenkins et alJ 12001 



[ Sheth et alJ200lHSeliak & Warrenl2004lTinker et alJ 2005b ; 
Warre net alJ I2005I) . In our approach, the parameters that 
determine the relationship between halo mass and the aver- 
age number of galaxies within a halo are constrained by ob- 
servational measurements of the projected two-point corre- 
lation functi on w„(r n ) (see, e.g., IZhengll2004t IZehavi et alJ 
120041 120051 lTinkeretalJl2005bl) . This technique is simi- 
lar to the condit ional luminosity function (CLF) analysis of 
Yang et al. (2003), but utilizes spatial clustering over a wide 
range of scales, linear to non-linear. Once the parameters 
of the occupation function have been determined, predictions 
can be made for other clustering statistics either analytically 
or through the creation of mock galaxy distributions with nu- 
merical simulations. In this paper we obtain the PVD by 
populating the halos of a high-resolution collisionless N-body 
simulation with galaxies such that both the spatial clustering 
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as a function of luminosity and the galaxy luminosity function 
match observations. This approach comp lements the analytic 
model for the PVD of lSlosar et al.l ( 120061) . and our qualitative 
conclusions about the PVD are consistent with their analysis. 

Directly measuring peculiar velocities requires requires 
knowledge of the true distance to galaxies, thereby remov- 
ing the Hubble flow from the measured redshift. This lim- 
itation subjects these measurements to shot noise and sam- 
ple variance. To circumvent this problem, many studies have 
used measurements of the two-point galaxy correlation func- 
tion in redshift space ^(r p ,r w ), which is dist orted along the 
line of sight by the pecu liar velocities (e.g., Davis & Peebles 
119831 lBeanet alJH98l . In this method, the galaxy PVD 
is obtained by modeling the redshift-space correlation func- 
tion at a given projected separation r p , assuming that the 
PVD is a constant along the line of sight. This method has 
be en used in more r ecent measurements of t he galaxy PVD 
bvlJing et alJ (fl^98l) . IZeha"vi et alJ J20021) . and lHawkins etail 
J2003I) . We will compa r e our predictions to the 2dFGRS 
results of lHawkins et al.l J2003I) . demonstrating qualitative 
agreement between the measurements and a model normal- 
ized to a matter density of fl,„ = 0.2. A quantitative com- 
parison is difficult d ue to the flux-limited sample used in 
lHawkins et aljj2003l) . and to possible systematic errors of the 
model used to infer the PVD from t he observed £,{r p , r„). 

Recent ly, Jing & Borner (2004, hereafter JB) and Li et. 
al. (2005, hereafter L05) have presented measurements of the 
PVD obtained by modeling observations of the redshift-space 
correlation function for volume-limited samples in the 2dF- 
GRS (JB) and Sloan Digital Sky Survey (L05). One notable 
result of these two studies is the luminosity dependence of the 
PVD at small scales. Both groups measure a clear trend of de- 
creasing velocity dispersion with increasing luminosity, until 
approximately M+- 1, whereupon the PVD rapidly rises. JB 
claim that this luminosity dependence is counterintuitive, and 
both papers show it is not reproduced in a halo occupation 
model. 

The technique of JB and L05 is to infer the PVD 
from redshift-space clustering in Fourier space using the 
"dispersion model" of the red s hift-space power spec- 
trum (e.g.. [Cole et alJll994l [i99l iPeacock & DoddslU994l 
lHatton & Coldl 1 99811 199lT see lScoccimarrol2004l for a thor- 
ough review and critique of this model). In their implemen- 
tation of the dispersion model, JB and L05 use the redshift- 
space power spectrum at a fixed £-scale and fit the data as a 
function of angle, as opposed to a constant projected separa- 
tion. The meaning of the velocity dispersion obtained in this 
case is less clear. JB argue that the qualitative transformation 
between the PVD in Fourier space and that in configuration 
space should be k «-> 1 /r, an inference motivated by the transi- 
tion between one-halo and two-halo clustering, which appears 
to occur at 1 Mpc" 1 h in £-space and 1 /r'Mpc in configuration 
space. The k «-> 1 jr transforma tion is also su pported by tests 
with mock galaxy distributions (Jing & Borner 2001). JB and 
L05 compare th eir results t o an implementation of a CLF fit to 
the 2dFGRS o flYang et al1d2003b . This CLF model predicts a 
PVD that increases with luminosity, nearly orthogonal to the 
observed trend. 

In this paper we test the JB implementation of the disper- 
sion model against our HOD models. We find that the fourier 
space dispersion measured by JB on k = 1 Mpc" 1 h can, in 
its interpretation as a configuration space dispersion on r = 1 
/i"'Mpc, be systematically off the true configuration spcae 
value at the 20-40% level, and only roughly follows the true 
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FIG. 1. — Panel (a): Satellite occupation functions for all nine magni- 
tude bins, as listed in Table 1 . The curves go in order of faintest to brightest 
galaxies from left to right. The thick gray line shows the total satellite occu- 
pation function for all galaxies. Panel (b): Central occupation functions for 
the same magnitude bins. The ordering of the lines is the same as in panel 
(a). The thick gray line is the sum of all central occupation functions. 



dependence of the PVD with configuration space scale. We 
pay particular attention to the question of what physical scale 
the k = 1 Mpc" 1 h measurements are probing. In configura- 
tion space, our HOD models predict that a v should be nearly 
independent of luminosity at r = 1 /r'Mpc. At smaller and 
larger scales, <j v has a weak dependence on luminosity, such 
that brighter galaxies have lower dispersions. Applying the 
dispersion model technique to our mock galaxy distributions, 
we demonstrate that the k=\/r comparison is reasonable for 
only a small set of luminosity samples, and in general there 
is no consistent comparison between dispersion model results 
and the true PVD at a given scale. Thus in Fourier space the 
dependence of the PVD on galaxy luminosity is steeper than 
that in configuration space, leading to the strong trend mea- 
sured in the JB and L05 results. 

2. DATA AND MODELING 
2.1. 2dFGRSData 

We make predictions for the PVD and its dependence on 
luminosity by inferring the mean halo occupation function 
from measurements of the projected galaxy correlation func- 
tion w p (r p ) of the 2dFG RS. These measur eme nts are s imilar 
to those presented in Nor berg et al.l ( 120011) and Norber g" st alJ 
(2002a), who measured w p (r p ) from a sample of ~ 160,000 
redshifts, but our data have been updated to include the full 
data release of the 2dFGRS, which contains ~ 221,000 red- 
shifts. Further details a bout the survey can be f ound in 
Colles s etalJ(|200lll200l and lNorberg etail (I2002bl) . 

The details of the clustering measurements can be found in 
Norberg et al. (in preparation). We present here a brief sum- 
mary of the calculations. Using the full 2dFGRS survey we 
create volume limited samples, with faint limits from Mb, = 
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-17.0 to Mbj = -21 .0, each sample 0.5 magnitudes 5 wide, for 
a total of nine samples. All galaxies brighte r than My, , = —21 
are gro uped into a single sample. As in Norberg et al. (2001, 
2002a), a careful account of the selection function is made 
and the correlation fun ctio ns are obta i ned us ing the standard 
lLandv & SzalavHl99 3) and Hamilton ( 1993) estimators, with 
typically 100 times more randoms than galaxies. The pro- 
jected correlation function is estimated by integrating £{r p , r n ) 
out to fV max = 70 /i _1 Mpc, providing a stable estimate for 
w p (r p ) out to at least r p = 40 /i _I Mpc. Due to the sensitiv- 
ity of the results on close pair incompleteness, clustering re- 
sults from 2dFGRS are only really reliable beyond typically 
~ 150 /r'kpc, which is the limit adopted in this paper too. 
The correlation function is measured in thirteen radial bins, 
spaced evenly by 0.2 in log 10 r beginning at log 10 r = -0.7. 

The errors on the clustering measurements are estimated by 
a bootstrap resampling technique on roughly equal sized sub- 
regions, of which there are 16 in total (8 in each 2dFGRS 
region, covering in each region approximatively the same sur- 
vey area). The main reason for choosing this technique over 
mock surveys is that it allows error estimates to be obtained 
for samples with different clustering properties, whereas cur- 
rent available mocks need to be tuned for each sample. More- 
over the existence of large structures in the survey (see, e.g., 
Baugh et al. 2004) are difficult to account for in the available 
mocks, but with bootstrap resampling of large subregions this 
type of cosmic variance can, to some extent, be account for. 
We create 100 bootstrap resamplings, which is sufficient to 
obt ain a stable error covarian ce matrix for each subsample 
(see Porciani & Norberg 2006 for details on stability tests for 
error matrices from bootstrap resamplings). We note though 
that this number of bootstrap resamplings is not large enough 
to address the covariance between the subsamples. For that 
reason, in this paper, we have to ignore the existing covariance 
between 'neighboring' volume limited samples. We use prin- 
cipal component analysis of the covariance matr ix to calcu- 
late th e x 2 °f a given HOD model (e.g. JPorciani & Giavaliscol 
2002), using the first twelve principal components. 

2.2. HOD Models 

To m odel w n (r n ), we use t he analytic mod el described in 
detail inlTinker e t alJ (l200 5bl) (see also. | Berlind & Weinberg ! 
l2002UCoorav & Shethll2002tlZheng|l2004l) . The mean occu- 
pation function is divided into two components; galaxies lo- 
cated at the center of mass of the host halo, (N ce n}M, and satel- 
lite galaxies distributed throughout the host halo with (usu- 
ally) the same distribution as the dark matter, (N S3 i)m- For 
galaxy samples defined by a luminosity or magnitude thresh- 
old, (N cen )M can be defined by a single halo mass scale below 
which halos are not able to contain galaxies in the sample. For 
magnitude bin samples, (N cen )M must have both a minimum 
and a maximum mass scale, the latter being defined by the 
minimum mass of the next brightest sample. A simple exam- 
ple would be to use a square window as the central occupation 
function, such that (N cen )' M = 1 for < M h < M^„, where 
; denotes magnitude bin. This parameterization assumes a 
unique mapping of central galaxy luminosity onto halo mass. 
A more physical model takes into accoun t the scatter between 
mass and luminosity, such as that used by Zheng et al. (2005) 
to model the occupation functions of semi-analytic and hydro- 

5 All magnitudes quoted in this paper assume h = Hq/ 100 = 1 . All magni- 
tude bins will be referred to by their lower limit or by a reference number i, 
as listed in Table 1 . 



dynamic simulations of galaxy formation. Here we modify 
their (N cen )M to take into account binned magnitude samples, 
i.e. 
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where (N ci:n )' M represents central galaxies in magnitude bin i, 
M' min is the cutoff mass scale for central galaxies, a\ M con- 
trols the width of the cutoff mass range, and erf is the error 
function. For luminosity threshold samples, M m \ n is defined 
as the mass at which (N ceR )M = 0.5, but in equation Q this 
mass can differ from M' min . The form of equation Q guar- 
antees that the sum of (N cen )' M over all i will never be larger 
than one, which would be unphysical. The form we use for 
the satellite galaxy occupation function is 
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where M' cut is the cutoff mass scale at which (N sM )' M transi- 
tions between a linear dependence of N on M and an expo- 
nential cutoff, M[ is the amplitude of the occupation function 
(forMj ^> M' cnV M\ is the mass at which halos host on average 
one satellite of magnitude i), and a sat is the power-law depen- 
dence of (N s - dl )' M on host mass. We assu me the scatter about 
the mean follows a Poisson distribution (Z heng et alJl2005t 
Kravtsov et al. 2004). Although a sat is a free parameter, we 
assume that it is constant for all magnitude bins. Allowing 
M' cut to vary as a function of magnitude allows enough flexi- 
bility in the satellite occupation function to model w p (r p ) for 
each bin. 

For spatial clustering, differentiating between central and 
satellite galaxies is important because pairs of galaxies that in- 
volve a central and a satellite galaxy have a different distribu- 
tion in the small-scale (i.e. 'one-halo') regime than satellite- 
satellite pairs. This differentiation is also important when 
modeling galaxy velocities. Central galaxies are assumed to 
be at rest with the center of mass of the halo, while satellites 
have random motions roughly equal to the virial dispersion 
of the dark matter halo. Thus, pairs of two satellite galax- 
ies in a common halo have a dispersion \/2 times larger than 
pairs of central and satellite galaxies. For two-halo pairs, in 
which each galaxy comes from a separate halo, the pairwise 
velocity of the halos contributes to all pairs, but the internal 
halo dispersions make an additional contribution to pars in- 
volving satellite galaxies. The ratio of satellites to total galax- 
ies is therefore important in determining the behavior of the 
PVD as a function of luminosity. Although the assumption 
that central galaxies have no random motions is a working 
approximation, in tests we find that the results of this assump- 
tion are nearly identical to models in which central galaxies 
h ave velocities in agreement w ith recent observational results 
of lvan den Bosch et alj j2005). who infer er cen sa 0.25<7 sat . For 
pairwise statistics these dispersions add in quadrature, thus 
central motions of this magnitude only increase the central- 
satellite dispersion by 3%. 

To calculate w p (r p ) for a given set of HOD parameters, we 
must first assume a cosmology and a form of the linear mat- 
ter power spectrum. Because we will use the results of our 
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fitting procedure to populate the halos identified in a colli- 
sionless N-body simulation, we take the same values as those 
used in the simulation. We assume a flat ACDM model with 
fi m = 0.3 and a linear amplitude of fluctuations of a% = 0.9, 
where erg is the variance of matter fluctuations on an 8 /r'Mpc 
scale. The shape of the l inear matter power spectru m was 
computed using cmbfast (Seliak & Zaldarriaaa 1996), with 
h = 0.7 and fib = 0.04. The simul ation was performed using 
the Hashed Oct-Tree code of Warren & Salmon ( 1993J). The 
simulation box is 400 /r'Mpc per side, with a total of 1280 3 
particles, giving a mass resolution of 2.5 x 10 9 ft Mq. Ha- 
los are identifie d in th e simulation by the friends-of-friends 
technique dDavis et alJ[l9 85) with a linking length 0.2 times 
the mean interparticle separation. For the lowest luminosity, 
Mi,j ~ -17 galaxies, the minimum mass roughly corresponds 
to a 50-particle halo. 

To calculate the model w p (r p ) for each bin, the HOD pa- 
rameters are set to match the abundance of galaxies within 
each magnitude bin n l „ as determ ined by the bj luminosity 
function of Norberg 'et al.l fc002bl) . with Schechter function 
parameters $* = 0.0161, M* = -19.7 and a = -1.21. This is 
done by calculating the value of M' min required to match n' 
once the other parameters for bin i have been chosen. For 
high-luminosity bins (i > 7, orM/^ < -20, which we will re- 
fer to as the 'bright' samples, while referring to i < 7 as the 
'faint' samples), we allow a[ M to be a free parameter. Recall 

that the shape of the cutoff in (N csn )' M physically represents 
the scatter between halo mass and central galaxy luminosity; 
at bright magnitudes this scatter is expected to be larger due 
to the steepness of the halo mass function in this regime and 
the errors in 2dFGRS galaxy magnitudes. For faint samples 
the scatter is expected to be smaller and magnitude errors do 
not increase it. We therefore fix cr[ M to be 0.15 for these 

samples. Leaving a[ M as a free parameter for these samples 
does not result in a statistically better fit to w p (r p ), and in tests 
we find that the results for the PVD are not sensitive to the 
value chosen. 

We find the best-fitting HOD models using the Markov 
Chain Monte Carlo (MCMC) technique (see, e.g., 
iDoran & MuelleJ 120041) . In MCMC, an HOD parameter 
set is randomly chosen from a Gaussian distribution centered 
on the last accepted parameter set in the MCMC chain. Each 
element in the chain contains the full set of HOD parameters 
for all nine magnitude bins. A model is automatically 
accepted into the chain if it results in a lower \ 2 value than 
the previously accepted model, and accepted with the prob- 
ability p = exp(-Ax 2 /2) if the \ 2 is higher. The number of 
parameters in our model is 22, consisting of two parameters 
(M\ and M' cut ) for each ; = 1-6, three parameters for each 
i = l— 9 (M'j, M' cat , and a[ gAf ), and a sa t, which is the same for 
all i. Because the values of M' min are set by the galaxy number 
densities, they are not free parameters 6 . The total \ 2 f° r eac h 
model is the sum of the \ 2 values obtained for the w p (r p ) in 
each magnitude bin. This large number of parameters makes 
MCMC an ideal technique for exploring the likelihood of 
possible solutions and of predictions these solutions make 
for the luminosity-dependence of the PVD. To check the 
chain for convergence we use the power spectrum method 

6 We neglect the errors on the galaxy number densities, which are of order 
a few percent. Due to the steepness of the halo mass function, these errors 
have little effect on M^ jn . The dominant source of uncertainty are the w p (r p ) 
measurements. 



of iDunklev et alJ (|2P05). The element with the minimum 
X 2 is taken to be our fiducial model, which is listed in Table 
1. Although we list the \ 2 value for each w p (r p ) in Table 
1, it should be noted that this model was chosen by its total 
X 2 - Therefore, for each individual bin i there will be models 
for which xj is lower than that listed in Table 1 . We select 
a subset of 1,000 models from the chain with which to test 
the constraints on HOD parameters and the PVD inferred 
from them. Table 1 also lists the galaxy bias factors b' of the 
fiducial simulation, as well as the dispersion in b' g from the 
1,000 MCMC models. The galaxy bias is straightforward to 
calculate analytically once the HOD parameters are known, 
and is determined by 

f 1 f°° dn , >; 
b' g =- / dM —— (N) M b(M), (3) 
s n' g J dM 

where dn / dM is the halo mass fun ction (iJenkins et alJ[2001l) 
and b(M) is the halo bias function dTinker et alJl2005bl) . The 
variance of b' g is ~ 1 -2% for most bins, with the tightest 
constraints for bins / = 4-7 (-18.5 > >-20.0). For lower 
luminosity bins the limited volume of the samples increases 
the errors on w p (r p ), while for higher luminosity bins the low 
number density of galaxies creates significant shot noise in 
w p (r p ). It should be stressed that the errors in V are for a 
given cosmological model. The galaxy bias obtained from 
w p {r p ) is degenerate with the assumed linear matter power 
spectrum P\m(k). 

Figure n shows the occupation functions from the fiducial 
model for all nine magnitude bins. Panel (a) plots, from left 
to right, (N sat )' M as a function of logM/, for i= 1-9. The thick 
grey line is the sum of all (Ny dt }' M , which effectively repre- 
sents the satellite occupation function for a magnitude thresh- 
old sample for galaxies brighter than Mbj = — 17. As expected 
from both the shape of the galaxy luminosity function and 
from the increased relative bias as a function of luminosity 
(see, e.g., Zehavi et al] l2005l [Norberg et al. 2002a), the am- 
plitude of the satellite occupation function M\ monotonically 
increases with the magnitude of the bin. The best-fit values 
of M' cut are not monotonic with luminosity, and for i = 1 and 
2 the satellite cutoff occurs at higher mass than i = 3. Due to 
the larger errors on w p (r p ) for lower luminosity bins, and the 
fact that (N SSL i)' M is not well-constrained for occupation num- 
bers significantly less than unity, there is a large range of M' cut 
allowed by the data for these magnitude bins. We will discuss 
this further in the following section. 

Figure^) plots (N cen )' M as a function of logM/, for the same 
magnitude bins. The thick grey line, once again showing the 
sum of all nine (N cen )' M functions, resembles the central occu- 
pation function of a magnitude threshold galaxy sample. For 
the brighter samples, where a{ M is allowed to vary, the shape 
of the cutoff becomes progressively softer with magnitude. 
This is expected for two reasons. First, the halo mass function 
contains an exponential cutoff above M/, ~ 10 13 h~ l M , and 
this steepness is expected to increase the physical scatter be- 
tween halo mass and central galaxy luminosity . HOD models 
of the brightest SDSS galaxies by Zehavi et al J 1120 05) that uti- 
lized an exponential cutoff were statistically preferred to mod- 
els with a step-function form of (N cen )M- Second, the magni- 
tude errors in the 2dFGRS are significantly larger than those 
of SDSS measurements. The exponential cutoff in the galaxy 
luminosity function creates asymmetric scatter, as galaxies of 
lower luminosity are scattered into a brighter bin more often 
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FIG. 2. — Panel (a): Model and observed correlation functions for the galaxies with -18 > M/,j > -18.5 ((' = 3 in Table 1). Points with error bars show 
the measured w p (r p ) from the full 2dFGRS data release. The solid line plots the best-fit analytic model, as listed in Table 1. The dotted line shows the w p (r p ) 
obtained from the numerical simulation, once the halos have been populated with the same HOD parameters. The dashed line shows the best-fit HOD when the 
parameters of Sanchez et al. 1 2005) are used to calculate the linear matter power spectrum. Panels (b) — (c) are the same as (a), but for magnitude bins i = 5, 7, 
and 9. 



than higher luminosity galaxies are scattered down. 

2.3. Mock Galaxy Distributions 

To create mock galaxy distributions, we populate the iden- 
tified halos in the N-body simulation. The number of satel- 
lite galaxies for each halo in each magnitude bin is chosen 
from a Poisson distribution with mean (N saJ )' M . The presence 
of a central galaxy is determined from a nearest integer dis- 
tribution from the sum of all (N cen )' M , and the luminosity of 
the central galaxy is chosen from the relative distribution of 
(N ci:n )' M . The central galaxies are placed at the center of mass 
of the halo, while the satellite galaxies are placed randomly 
th roughout the halo acc ording to the universal density profile 
of lNavarro et alJ ill 9971). For each halo mass, the appropriate 
parameters of the Navarr o et al.l (119971) profile are calculated 
for our assumed cosmology from the method of Bulloc ket al.l 
teOOll) . The number density profiles of the satellite galaxies 
are assumed to follow the density profile of the dark matter. 
As stated above, the central galaxy is given the velocity of 
the center of mass of the halo, while the satellite galaxies are 



given an additional random velocity in each direction drawn 
from a Gaussian distribution of width 



GM h 
2R 



(4) 
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where /?200 is the radius at which the mean interior density of 
the halo is 200 times the universal mean, a value that corre- 
lates well on average with the mean densities of friends-of- 
friends halos with a 0.2 linking length. The relation between 

M h and R 2 qo, R 2 m = (3M/,/47r200/9 m ) 1/3 , makes a h oc AfJ /3 '. 
Equation 10} assumes no velocity bias for satellite galaxies. 
A more general form of equation (|4) would have a constant 
of prop ortionality az relating satellite motions to dark matter 
(Berlind & Weinber gl2002HTinker et all2005al) . 

Figure |2l shows the w p (r p ) data and the HOD models for 
four of the magnitude bins, as labeled in each panel. The 
points with error bars are the 2dFGRS measurements, the 
solid lines are the fiducial analytic model, and the dotted lines 
show the measured correlation functions from the populated 
N-body simulation. The total \ 2 of this model is 105.3 for 
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TABLE 1 

HOD Parameters for Models in FigureQ 



; 




9 

x 2 


M miB 


M max 
'"ceil 


M\ 


M cut 




h 


K 


1 


-17.0 


15.6 


1. 11x10" 


1.27 x 10 11 


1.01 xlO 13 


4.17xl0 12 


0.15 


0.90 


0.02 


2 


-17.5 


6.5 


1.57x10" 


2.01 x 10" 


1.08xl0 13 


1.07X10 13 


0.15 


0.92 


0.03 


3 


-18.0 


6.3 


2.38x10" 


2.86 x 10" 


1.44xl0 13 


5.71x10" 


0.15 


0.94 


0.015 


4 


-18.5 


10.9 


3.71x10" 


5.12 x 10" 


1.80xl0 13 


5.64xl0 12 


0.15 


0.93 


0.011 


5 


-19.0 


13.1 


6.73x10" 


9.15 x 10" 


2.61 xlO 13 


1.84xl0 12 


0.15 


0.96 


0.011 


6 


-19.5 


12.5 


1.40X10 12 


2.32 x 10 12 


3.83X10 13 


1.45xl0 12 


0.15 


1.03 


0.007 


7 


-20.0 


9.7 


3.37X10 12 


4.66 x 10 12 


6.93X10 13 


6.49xl0 12 


0.09 


1.13 


0.009 


8 


-20.5 


25.5 


1.22x10" 


2.67 x 10 13 


1.82X10 14 


2.19xl0 13 


0.33 


1.29 


0.03 


9 


-21.0 


5.2 


6.05xl0 13 




2.49xl0 14 


3.48xl0 14 


0.56 


1.69 


0.08 



NOTE. — All masses are in units h 1 Mq. The best-fit value of o sa t from equation |2j is 1.03. M™J X is the mass at which {N C en)' M is 
maximum. b g is the galaxy bias factor, and Sj, is the dispersion of bias values from the MCMC chain. 



9 x 12-22 = 86 degrees of freedom, yielding a \ 2 P er de- 
gree of freedom of 1.22. But from inspection of Figure|2]it is 
apparent that the cosmology chosen lacks large-scale power 
relative to the data. If these data are fit using a Pn n (k) consis- 
tent with the recent analysis of the 2dFGRS l arge-scale g alaxy 
power spectrum and CMB results (Sanchez et al. 2005), the 
difference at large scales is ameliorated, and the total \ 2 is 
reduced. These fits are shown with the dashed lines in Figure 
|2] The fits at small scales are unchanged, but at large scales 
the models match the amplitude of the observations. The \ 2 
per degree of freedom for these fits is 0.61, a value that sug- 
gests the possibility that the errors on the data have been over- 
estimated, perhaps by neglecting the covariance in w p (r p ) be- 
tween luminosity bins, or by using too many principa l compo- 
nents in the error analysis (Porciani & Norberg 2006). We use 
the standard ACDM fits i n order to populate th e simulation. 
Using a simulation with a Sanchez et al. (2005) power spec- 
trum would likely produce small differences in the predicted 
PVD, but it would not radically change our predictions be- 
cause the best-fit halo occupation functions are not very sen- 
sitive to the assumed power spectrum shape. 

2.4. Satellite Fraction of Galaxies 

A key quantity for interpreting the results of the HOD anal- 
ysis is the satellite fraction / sat , determined for each bin from 
the HOD parameters by 



1 f°° dn , w . 



(5) 



We plot this quantity as well as the different mass scales from 
the occupation functions in Figure[5] Panel (a) plots the radii 
of M' cut halos as a function of magnitude for both the fiducial 
model and the MCMC chain. We present this quantity as a 
radius rather than a mass because the physical sizes of ha- 
los influence the scales at which one-halo and two-halo pairs 
can contribute. The solid black line represents R(M' cul ) for the 
fiducial model. The shaded region represents the inner 68% 
of the distribution of R(M' cut ) about the median value from the 
MCMC chain. The distribution of R(M' cut ) is calculated for 
each bin individually, so the upper and lower bounds on the 
shaded region do not represent a single element in the chain, 
but the shaded region can be taken as an estimate of the error 
in M' cut . The bin-to-bin variation seen in the fiducial model is 
common for individual models from the chain. The midpoint 
of the shaded area is not defined by a single model, thus the 
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FIG. 3. — Panel (a): Radii of M' cM halos, below which the satellite occu- 
pation function is exponentially suppressed, as a function of magnitude. The 
thick solid line is the relation from the fiducial model, while the shaded region 
represents the inner 68% of the distribution of M' cut values about the median, 
taken from the MCMC chain. Panel (b): Radii of M' mm halos (bottom shaded 
region), the cutoff scale for central galaxies, and M' ons , halos (upper shaded 
region), which on average contain one satellite. These results are taken from 
the MCMC analysis only. Panel (c): Fraction of galaxies that are satellites as 
a function of magnitude. Line and shaded region are as in panel (a). 



best-fit model need not follow it. The constraints on M' cal are 
weakest at lowest luminosities, where the observational errors 
are largest. The shaded region narrows for Mb, < -19.5 due 
to the smaller errors on w p (r p ). 

Figure |5J) shows the results for M' min and M' one for the 
MCMC models only, where M' olK is the mass at which 
(Nsat)'M = 1> which can be larger than M\ due to the expo- 
nential cutoff in equation 0. The shaded regions show the 
same distribution about the median as in panel (a), with the 
lower shaded region representing the results for M' min and the 
upper region representing M' one . The constraints on both these 
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FIG. 4. — Pairwise velocity dispersion as a function of radial separation 
for four magnitude bins, measured from the populated simulation. 



parameters are significantly tighter than those on M' cuV The 
majority of one-halo pairs come from halos with ~ 1 satellite 
galaxy, so w p (r p ) provides the best constraints on the occupa- 
tion function around M' one . The value of M' cut influences the 
shape of (N s - dt )' M where satellite occupation is typically less 
than ~ 0.5 and it is therefore less well constrained by w p (r p ). 
The value of M' min is tightly constrained by both w p (r p ) and 
the galaxy number density. Not plotted is the distribution of 
a sat , which is nearly Gaussian, with a mean of (a sat ) = 1 .03 
and a dispersion of 0.03. It should be noted that this value 
of a sat is degenerate with the assumed value of <7g; a lower 
value of <jg would increase the inferred value of a sat in order 
to match the observed clustering. 

Figure |3j; plots / s ' at as a function of magnitude for both the 
fiducial model (solid line) and the MCMC chain (shaded re- 
gion). As with M' cuV the constraints on / s ! at are weaker for 
fainter samples, where it is constrained to ~ 25% about a me- 
dian value of approximately 0.24. The median value of / s ! at 
drops gradually to ~ 0.20 for ~ M* galaxies, then falls 
rapidly for samples brighter than -20, reaching / s ' at = 0.06 ± 
0.02 for i = 9. These trends of the satellit e fraction are in 
good agreement with the results of Mandelbaum et al. I J20051) . 
who derive satel lite fractions from meas urements of galaxy- 
galaxy lensing. Mandelbaum et al. (2005) report / sat « 30% 
forL <L* galaxies (in Sloan r-band) and/ sat ss 0.13 for galax- 
ies significantly brighter than L„. 

This estimate of / sat depends only weakly on the value of 
erg assumed. To match the observed w p (r p ) for a lower value 
of erg, a sat increases to compensate for the lack of high-mass 
halos. This increases the mean bias of galaxies, but for / sat the 
effects of higher a sat and fewer massive halos nearly cancel. 
In tests where we assume erg = 0.7, the median / s ' at is lower by 
~ 10% than the results in Figure [3] but the 68% confidence 
regions overlap significantly at all magnitudes. The shape of 
/ sa t' particularly the rapid turnover at M*, is also preserved. 



3. LUMINOSITY AND SCALE DEPENDENCE OF THE 
PAIRWISE VELOCITY DISPERSION 

3.1. Theoretical Predictions 

We calculate the PVD from our populated simulation by 
binning pairs in real space by their radial separation and cal- 
culating the line-of-sight pairwise dispersion as a function of 
r. We focus on a v as a function of three-space separation r 
rather than projected separation r p for two reasons: present- 
ing a v {r) offers a cleaner picture of the PVD at small scales, 
and the measurements of JB and L05 are intended to be of 
the PVD at fixed scale and not as function of projected sep- 
aration. Figure 0] shows the radial dependence of <j v for four 
bins of increasing magnitude (samples i = 2,4,6,8 in Table 
1). The general sha pe of the curves is sim ilar to those found 
for dark matter (e.g., Efstathiou et a l.l!988l) : at large scales, <r v 
has little scale dependence, but at r ~ 3 /r'Mpc, a v rapidly in- 
creases towards smaller r, peaking at r ~ 1 /r'Mpc. This fea- 
ture represents the transition between two-halo pairs, which 
are dominated by pairs of central galaxies, and one-halo pairs 
dominated by pairs of satellite galaxies. At r < 1 /r'Mpc, the 
PVD turns over and monotonically falls with smaller separa- 
tion. This feature reflects two aspects of the one-halo term. At 
smaller separation, the relative contribution of central-satellite 
pairs to the overall pair count increases. Furthermore, as r de- 
creases, lower-mass halos begin to contribute one-halo pairs, 
which further reduces the PVD. The former effect is more pro- 
nounced for brighter samples because of the smaller fraction 
of satellite galaxies for these samples. 

At a given separation, the trend of a v with luminosity can 
be inferred from Figure 0] For r > 1 /r'Mpc, the dependence 
of the PVD on magnitude varies from roughly flat to a mod- 
erate decrease with increasing luminosity. At smaller separa- 
tions, the dependence with luminosity becomes stronger. The 
location of the transition between between one-halo and two- 
halo pairs also depends on luminosity, with lower-luminosity 
galaxies transitioning to one-halo pairs at larger separation 
(~ 3 /r'Mpc for i = 2) as compared to bright galaxies (~ 1 .5 
/r'Mpc for i = 8). These results are in good agreement 
with th e analytic HOD calculation of the PVD in Slos ar et all 
(2006); both models predict the same characteristic shape to 
CT v (r), as well as the same behavior with luminosity at r < 0.5 
/r'Mpc. There are minor differences between the predictions, 
such as the trend of the PVD with luminosity at r > 1 /r'Mpc, 
which in their model increases with increasing luminosity, 
which is the opposite of the trend seen in Figure|4] This dis- 
crepancy results from Slosar et al. (2006) fixing the satellite 
fraction to be a constant value of 0.2 with lumin osity. Here 
we ob tain / sat directly from the clustering data. Slos ar et all 
(2006) conclude that the amplitude of the PVD is highly sen- 
sitive to the value of / sat , and we concur. 

We use the results of our fiducial model to interpret the lu- 
minosity dependences shown in the previous figure. Figure|5] 
breaks down the relative contributions of the different types 
of pairs as a function of magnitude within the radial bins at 
0.1, 1, and 3 /r'Mpc. Within isolated halos, galaxy pairs can 
be either central-satellite and satellite-satellite. These types of 
pairs can also exist between two distinct halos, in addition to 
the central-central contribution. For brevity, we refer to these 
pair types as cen-sat, sat-sat, and cen-cen in the Figure leg- 
end and in the discussion below. Panel (a) plots a v for each 
type of pair at r = 0.1 /r'Mpc, as well as the overall PVD 
with the solid line. At this separation, all galaxy pairs are 
one-halo. The error bars on the solid line are obtained from 
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jackknife sampling of the simulation volume into eight oc- 
tants. Each octant is roughly equal to the volume of the M* 
sample used here. The shaded region once again shows the in- 
ner 68% about the median value of a v at each luminosity from 
the MCMC chain. For faint galaxies, a v is entirely accounted 
for by the sat-sat contribution because (N cen )' M and (iVsat)jtf f° r 
these samples do not overlap in mass; halos massive enough to 
host satellites are too massive to have a low luminosity central 
galaxy. At magnitudes brighter than M[ n = -19, a v decreases 
significantly even though the sat-sat PVD remains constant. 
For the brightest sample the overall PVD is equal to the sat- 
een dispersion. Panel (b), which plots the relative contribu- 
tion of each pair type to the overall pair count as a function of 
magnitude, shows why this occurs. Below Mtj = -19, sat-sat 
pairs account for all pairs. At ~ M*, the contribution of cen- 
sat pairs rapidly increases due both to the presence of central 
galaxies and to the lower fraction of satellites. For M/ v < — 21, 
all pairs involve a central galaxy. For this sample, (N sat )' M = 1 
at 4 x 10 14 /r'M©. Although the abundance of halos more 
massive than this value is low, they do exist in the simula- 
tion volume, making it possible to have sat-sat pairs at this 
magnitude. But because the volume of these halos is rela- 
tively large, the probability of having two very bright satellite 
galaxies within 0.1 /r'Mpc is low. 

Figure^ plots o> for each type of pair at r = 1 /i~'Mpc. At 
this separation, which combines both one-halo and two-halo 
pairs, there are contributions to a v from all three pair types. 
However, the overall PVD lies between the sat-sat and cen- 
sat dispersions, indicating that one-halo statistics still domi- 
nate. This is confirmed in panel (d), which plots the relative 
fraction of the pair counts. Sat-sat pairs dominate at all mag- 
nitudes until the brightest bin. For Mb, < -21 galaxies, the 
contribution of cen-sat pairs is equal to sat-sat pairs. Fig- 
ure [5Ji showed that the radius of M' cut halos for this sample 
is > 1 /i~'Mpc, roughly equal to the radial separation shown 
here. For faint galaxies, cen-sat pairs are more common than 
pairs between central galaxies, but the abundance of cen-sat 
pairs gradually decreases with magnitude and cen-cen pairs 
become more common at M*, a result of the lower / sat val- 
ues. Panel (e) plots a v for each pair type at r = 3 /i _I Mpc. 
At this separation two-halo pairs dominate the statistics, and 
the overall PVD lies between the cen-sat and cen-cen values. 
In panel (f), cen-sat and cen-cen pairs are nearly equal for 
faint galaxies, but there is once again a transition at M* where 
the contribution of cen-cen pairs increases and the other pair 
types become less common due to the lower satellite fraction. 

At small separation, well within the one-halo term, the de- 
pendence of the PVD on galaxy luminosity is strong, with 
brighter galaxies having lower dispersions (with the excep- 
tion of the brightest galaxies, which live in the highest-mass 
halos). At intermediate separation, where the transition be- 
tween one-halo and two-halo pairs occurs, the dependence on 
luminosity is washed out due to the dominance of satellite 
pairs at all magnitudes. Outside the one-halo term, but still 
in the non-linear regime, there is a weak dependence of a v 
on luminosity, falling from ~ 700 kms" 1 for Mb, = -17 to 
~ 500 km s" 1 for Mbj = -20.5. (From Figure0]it can be seen 
that this slope becomes shallower at even larger scales.) At 
all three scales, the results of the MCMC chain demonstrate 
that these predictions are robust given the errors in the w p (r p ) 
measurements. For the latter two separations, the width of the 
shaded region is comparable to the jackknife errors from the 
simulation itself. At face value, these results are at odds with 



the observational results of L05 and JB, which show a strong 
dependence of a v on luminosity at k = 1 Mpc" 1 h, which they 
argue is equivalent to r = 1 /r'Mpc. Our HOD results are also 
inconsistent with the CLF prediction shown in both papers, 
which shows a strong trend of increasing PVD with luminos- 
ity. In section §4 we will use the JB technique of obtaining 
the PVD from the redshift-space power spectrum to make a 
proper comparison of our HOD predictions with these results. 

3.2. Comparison with 2dFGRS Measurements 

Estimates of the PVD from the redshift-space correlation 
function do not show the shard transiti on at r ~ 2 /r'Mpc 
seen in FigureEl Results fromlJing et alJ (119981) . IZehavi et all 
(2002), and Hawkins et al. (2003), from the Las Campanas 
Redshift Survey, the SDSS, and the 2dFGRS, respectively, 
show a smooth increase in PVD from large to small scales. 
However, as discussed in § 1 , the PVD is not measured directly 
from the data but inferred by model fitting, and the quantity 
that is fitted can be interpreted as an "effective" velocity dis- 
persion along the line of sight, calculated as a function of pro- 
jected separation r p . In Figure [6] we plot the mean line of 
sight <t v as a function of r p for pairs with a real-space line 
of sight separation of less than 40 /r'Mpc. The four curves 
are for the same magnitude bins as shown in Figure [4] For 
comparison, the data from Figure |4] for i = 2 are shown with 
the open squares. The sharp transition from one-halo to two- 
halo pairs is smoothed out when plotting a v (r p ) because at a 
given r p pairs with larger separations (and lower dispersion) 
contribute to the average. 

Also shown i n Figu re [6] are the PVD measurements from 
lHawkins et al.l d2003l) . These data follow the axis on the 
right-hand side of the plot, the values of which are lower 
by a factor of 0.8. Thus the amplitude of our HOD pre- 
diction is higher than the data. These calculations assume 
ft,,, = 0.3 and no velocity bias. Both of these assumptions 
could be wrong, and there exists growi ng observational ev- 
idence indicating a lower v alue of (Bahcall et alJl2003l 
I van den Bosch et alJ I20031 ICole etail I2005I [Sanchez et alJ 
I2005I lMohavaee&Tullvll2005l iTinker et alJ l2005b!) and 
some theo retical evidence for moderate negative velocity 
hiasdKlvnin etal.ll 1999: Berlin7^l2003l fYoshikawTeTatl 
2003; Slosaretal. 2006). Therefore we consider the over- 
all normalization of our prediction to be somewhat arbitrary. 
Read against the right-hand axis, the curves represent mod- 
els scaled 7 by a factor of (0.2/0.3)° 55 = 0.8 to depict a uni- 
verse with fi m = 0.2. For fixed power spectrum shape, chang- 
ing Sl m does not a lter the shape of the halo mass function 
( Zhen get alJl2002l) . Thus the occupation functions in Fig- 
ure once scaled to the new value of Q, m , produce identical 
fits to w p (r p ) as in Figure[2] Once scaled to this new value of 
f2 m , the amplitude and scale dependence of the PVD from the 
HOD models is generally consistent with the 2dFGRS data. 
There is some discrepancy at r p = 0.1 /r'Mpc, but this data 
point suffers from fiber collision effects and incompleteness 
due to blending of objects in the input catalog. A quantita- 
tive co mparison between the data and the models is difficult 
because Hawki ns et al.l (12003 ) analyze a flux-limited sample, 
and we have not replicated the details of their modeling pro- 
cedure. However, a flux-limited "prediction" would be a lin- 
ear combination of the different model curves (and the cross- 

7 Halo pairwise velocities scale as CjJ, 6 even in the non-linear regime 
1 Zheng et al. 2002), while virial velocities scale as Hj], 5 , thus to a good ap- 
proximation the pairwise velocities scale as fij}; 55 at all scales. 
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FIG. 5. — Panel (a): PVD at r = 0.1 /i Mpc, broken down by pair type, as a function of luminosity. The dashed line is the PVD for satellite-satellite pairs 
only, the dash-dotted line represents central-satellite pairs, and the dotted line (not seen in this panel) is for central-central pairs. The solid line is the PVD for all 
pairs. All curves are obtained from the Fiducial HOD model in Table 1. The shaded region is the inner 68% of the distribution of a v values about the median 
in each magnitude bin, obtained from the MCMC chain. Panel (b): The fraction of the total pair count contributed by each type of pair. Line types are the same 
as in panel (a). Panels (c) and (d): PVD and pair fraction for each pair type at r = 1 /r'Mpc. Line types are as above, with the dotted line representing the 
central-central contribution. Panel (e) and (f): Same as above, but for r = 3 /r'Mpc. 



correlated velocity statistics). All the (scaled) model curves 
show a smooth increase in a., (r p ) from - 400 km s _I to - 600 
kms" 1 , in qualitative agreement with the data. 

4. INTERPRETING THE RESULTS OF THE DISPERSION 
MODEL 

JB and L05 measure velocity statistics for volume-limited 
samples, enabling easier comparison to theoretical predic- 
tions. However, the quantity that they measure in not the PVD 
itself but the parameter at of the dispersion model, which con- 
volves linear theory with an exponential distribution of ran- 
dom galaxy velocities. With this model, the redshift-space 
power spectrum P z {k,jj) is related to the real-space power 
spectrum //) by 

P z {k^) = P R {k) (1 +/3 M 2 ) 2 (1 + alii 2 k 2 /2)~ l , (6) 

where /i is the cosine of the angle of the wavevector with the 
line of sight, (3 = il^ n 6 /b g , and er* is the non-linear galaxy ve- 
locity dispersion. The exponential distribution assumed in 



equation O is both observa t ionally and theoretically moti- 
vated JDavis & Peeblesll983tlShethll9 96). In a recent paper, 
Scoccimarro (2004) presents a thorough critique of the dis- 
persion model, pointing out several deficiencies. First, the 
physical probability distribution of random velocities implied 
by equation (|6j is not a true exponential dispersion, but in 
fact contains a discontinuity at the mode. Second, the veloc- 
ity distribution function of dark matter deviates significantly 
from an exponential at some scales. T his is also true of halo 
pairwise velocities ( Zure ket all 1 9941) . Inside the scale of the 
one-halo term, an exponential is a good approximation for 
galaxies, but in the two-halo term some deviations are to be 
expected. Third, equation (|6j is also only a true convolution 
if the parameter is assumed to be constant with scale, but 
this assumption is contradicted by the results of Figure|4] The 
scale dependence of the PVD is part what JB and L05 wish to 
probe, thereby removing one of the theoretical underpinnings 
of equation In light of this, it is necessary to quantify the 
accuracy of the JB implementation of the dispersion model in 
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FIG. 6. — Pairwise velocity dispersion as a function of projected separation for the same four magnitude bins shown in the Figure|4](as a function of 3-d 
separation). Only pairs with a line-of-sight separation less than 40 /r'Mpc are counted. The open squares are the results from Figurel4lfor the faintest bin shown 
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universe with O m = 0.2 (see text for details). The filled circles show the 2dFGRS measurements for a flux-limited sample from Hawkins et al. 12003), which 
should be viewed with respect to the right-hand axis. 



order to assess the possible sy stematic errors accrued through 
its use. Uing & Bornerl (120011) test the dispersion model with 
mock galaxy distributions, finding favorable results, but they 
apply these tests to multiple realizations of one distribution 
and not as a function of luminosity. Here we test the disper- 
sion model against all nine magnitude bins. 

We calculate the correlation function in redshift space 
£(r p ,r„), where r p is projected separation and is line- 
of-sight separation. We measure ^{r p ^r^) for < r p < 50 
/i _1 Mpc, and -40 < r v < +40 /r'Mpc. We Fourier transform 
into Pz(k, y) by 

Pzikp,^) = 2tt J dr w J ^<ir ; ^(r p ,r Tr )r p cos(^ Tr r 7r )7 (fep?- p )Wj(s) 

(7) 

where k p and k n are the projected and line-of-sight directions 
in £-space, Jq(x) is the zeroth-order Bessel function, and 

W g (s) = ex V (?pj, (8) 

is a weighting function utilized to reduce the influence of large 
scale distortions, where s 2 = r^ + r 2 , and S = 20 /r'Mpc (see 
equations [6] and [8] in JB). We calculate the error in Pz(k, /i) 
from eight jackknife samples and fit equation at each k, 
determining the best-fit values of Pz(k, 0) and by x 2 mini- 
mization. We fix the value of j3 for each magnitude bin using 
the value of b g from Table 1, but in practice we find the re- 
sults are insensitive to the value of (3 (in agreement with JB 



and L05). 

Figure presents the results of the JB method as applied 
to our fiducial model. The filled circles plot the dispersion 
model results as a function of k, the solid lines plot the true 
PVD in configuration space (from Figure |4j as a function of 
ir/r, and the dotted lines plot the same curves but as a func- 
tion of 1 jr. The errors on the solid points are the dispersion 
between the three projections of the simulation volume. For 
all magnitudes, the dispersion model results show the scale 
dependence expected from Figure|4j at large scales (small k), 
there is an asymptotic value of crj, and at smaller scales the 
dispersion recovered from the equation rises, reaching a 
peak between k = 1 — 3 Mpc" 1 h for most bins. For i = 4 and 
5, ait has no clear maximum value and continues to increase 
monotonically with k. Qualitatively, the dispersion model cor- 
rectly reproduces the scale-dependence of the pairwise veloc- 
ity dispersion that is measured in configuration space. In de- 
tail, however, the comparison between o> and er,, is problem- 
atic. For the lowest three magnitude bins, the transition from 
two-halo to one-halo pairs is complete by k ~ 1 Mpc" 1 h, mak- 
ing a favora ble comparison with t he a v {\/r) (dotted) curves, 
as shown in Jing & Bornej fcOOll) . For brighter samples, this 
transition moves to larger k, such that for i = 6-8, <Jk clearly 
follows the a v (ir/r) (solid) curves. For the other samples, it is 
not clear that either line best represents the dispersion model 
results. The results of Figure underscore the difficulty in 
translating the results of the dispersion model at a given A: to a 
single physical scale. For faint galaxies, the asymptotic value 
of (T/c at small k is larger than the high-separation value of 
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FIG. 7. — Dispersion model results for the fiducial HOD model. Each panel represents a magnitude bin as labeled in top-left corner. Filled circles plot 
as a function of k, and solid lines represent the true PVD measured directly from the simulation (see Figure|4j plotted as a function of rc/r. Error bars are the 
dispersion about the mean of the three projections of the simulation volume. The dotted lines plot the same results as the solid lines, but as a function of 1/r. 
Note that the range on the y-axis has been increased to a maximum of 1,500 kms~' for the bottom right panel (i = 9), which plots the results for the brightest 
magnitude bin. 



a v by ~ 200 kms" 1 , or 40%. For samples near M„, a* and 
<7 V roughly agree at small k. For the brightest samples, is 
smaller than cr v at large scales. 

Figure [8] compares the measurements of JB and L05 at 
k = 1 Mpc" 1 h with the dispersion model results of the fidu- 
cial model from Figure The open squares and the solid 
circles are the JB and L05 results, respectively. The results of 
both JB and L05 show a clear trend of lower velocity disper- 
sion for increasing galaxy luminosity in the range M» + 1 > 
Mi,j > M» - 1 . The claim of JB is that this result is counter- 
intuitive; brighter galaxies should reside in higher-mass halos 
and therefore should be subject to larger virial motions. This 
claim is reinforced by the predictions of a conditiona l lumi - 
nosity function model of the 2dFGRS by Yang et al. (2003). 
The dotted line in Figure |8]plots the predictions of their CLF 
model, obtained from inspection of Figure 9 in JB and Fig- 
ure 7 in L05. The CLF prediction is nearly orthogonal to the 
observational results, with a monotonic trend of increasing a> 
with luminosity. 

The large solid squares plot the HOD dispersion model re- 
sults taken from Figure [7] For several of the magnitude bins, 
the value of <Jk at k = 1 is slightly off the trend shown by the 
other points, so we average the model results at k = 1 with the 
data points on either side. The error bars represent the error in 
the mean of the nine values (three projections for each of three 



values of k). In configuration space, there is little to no trend 
of Uy with luminosity at r = 1 /! _1 Mpc (Figure^). The disper- 
sion model results of the fiducial model show a strong trend 
with luminosity that is consistent with the dependencies seen 
in the JB and L05 measurements. Figure[7]demonstrates that 
at k = 1 Mpc" 1 h, cr/i for faint galaxies probes the regime where 
satellite pairs dominate the statistic, while for brighter galax- 
ies the k = 1 Mpc" 1 h scale is outside the one-halo regime, 
where the dispersion is lower. This transition produces an en- 
hanced luminosity dependence in au at this scale. The slope 
of the HOD results does not appreciably change if we plot 
<7jt at k = 1 Mpc" 1 h without averaging over the neighboring k 
values. 

As in Figure |5J the HOD predictions for this cosmological 
model are higher than the measurements. The solid line rep- 
resents the HOD model, scaled again by (0.2/0.3)° 55 = 0.8. 
This scaled version of our HOD prediction is consistent with 
the data in both slope and amplitude, even reproducing the 
sharp upturn of the PVD for the brightest galaxy sample. 
Jing & Borner (2001) applied the dispersion model to the Las 
Campanas Redshift Survey, also finding that a model with 
fl m = 0.2 was a better fit to the data than f2 m = 0.3 (for erg = 1). 

The expectation of JB that brighter galaxies should occupy 
more massive halos is borne out by our HOD modeling. But 
this creates a strong correlation between the luminosity of 
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central galaxies and halo mass, while for satellite galaxies 
the lowest luminosities dominate the distribution (by number) 
within a given halo. The PVD of faint galaxies is high be- 
cause of many of them are satellites in high mass halos, and 
the PVD of luminous galaxies is lower than naively expected 
because many of them are central galaxies of massive halos 
that do not inherit the halo velocity dispersion. The results 
of JB and L05 are thus naturally explainable in the halo oc- 
cupation context through the effect of / sat , discussed in the 
previous section, and by the systematic effects of using the 
dispersion model to estimate the PVD. 

Why are the predictions of two similar methods (e.g., the 
HOD and the CLF) so clearly at odds? The real-space clus- 
tering results of the JB implementation of the CLF overes- 
timate the observed power spectrum for k > 1 Mpc -1 h by 
a factor of 3-5 for galaxies brighter than M/, y = -19.5 (see 
their Figure 6). This indicates that the number of bright satel- 
lite galaxies is also overestimated, leading to more one-halo 
pairs at these magnitudes than is observed, as well as a higher 
satellite fraction at these magnitudes. A model that predicts 
too many bright satellite s would predi ct an increasing PVD 
with luminosity, as ISlosar et alJ (|2006) demonstrate. In the 
lYang et all d2003l) CLF, galaxy clustering data are used only 
at the scale length ro of the correlation function, well outside 
the one-halo term. Therefore, models that overpredict small- 
scale clustering are not disfavored by high \ 2 values in their 
analysis. JB make a s imilar point about the best-fit parame- 
ters of the Yang et all J2003I) model, noting that a model with 
a steeper faint-end slope of the CLF could be more consis- 
tent with the PVD measurements. The differences between 
our HOD results and the CLF come as some surprise because 
of the good agreement between the two methods in previ - 
ous analyses llTinker et alJ2005btlvan den Bosch et al.l2 003). 
These differences are therefore specific to the statistic be- 
ing modeled here. The utilization of small-scale clustering 
for determining the occupation functions makes our approach 
more robust for predicting other statistics at small scales. For 
statistics that require knowledge of the mean properties of the 
galaxies within halos, such as the mass-to-light ratio or the 
overall PVD of galaxies, there should be less bias in CLF re- 
sults. 

5. SUMMARY 

In this paper we use the halo occupation distribution frame- 
work to make robust predictions of the pairwise velocity dis- 
persion of galaxies through modeling of the galaxy two-point 
correlation function. Our technique is to constrain the param- 
eters of the halo occupation function (N) m by matching ob- 
servations of the luminosity-dependent projected correlation 
function w p (r p ), then populate the halos of an N-body simu- 
lation with mock galaxies drawn from the inferred HOD. In 
our analysis, we make a distinction between central galaxies, 
which we assume move with the center of mass of the host 
halo, and satellite galaxies, which we assume have a veloc- 
ity dispersion within each host halo equal to that of the dark 
matter. This distinction is crucial to our predictions because 
the pairs that involve central galaxies have a lower dispersion, 
so the fraction of satellites strongly influences both the lumi- 
nosity and scale dependence of the PVD in our predictions. 
At large separations, the PVD is largely scale independent 
because halo velocities vary weakly over the scales probed 
here. At r ss 2 /r'Mpc, the PVD rapidly increases as satellite- 
satellite pairs from massive halos dominate. At r < 1 /r'Mpc, 
the PVD decreases with smaller separation because central- 



satellite pairs become more common. 

We have focused on the luminosity dependence of the PVD 
in light of recent claims that the measured luminosity depen- 
dence is in conflict with the expectations from halo occupation 
models. JB and L05 measure a PVD that decreases with lumi- 
nosity, which is inconsistent with the expectation that brighter 
galaxies should have higher dispersions because they occupy 
only high-mass halos. We have shown that this expectation 
is incorrect; bright galaxies occupy the centers of high-mass 
halos, so their dispersions are not enhanced by the large virial 
motions within those halos. The dependence of a v on M/,, is 
strongly influenced by the fraction of galaxies that are satel- 
lites at each magnitude, a quantity that is ~ 25% for faint 
galaxies and rapidly declines for < M* . In the CLF pre- 
diction for the PVD shown in JB and L05, the PVD mono- 
tonically increases with luminosity. This prediction does in- 
clude the distinction between central and satellite galaxies, 
but it does not utilize w p (r p ) data at small scales to constrain 
the occupation function, so it cannot accurately constrain / sat . 
While we find that observational uncertainties in w p (r p ) and 
degeneracies among HOD parameters allow some range in 
our predictions of the PVD, none of our 1,000 MCMC mod- 
els predict a monotonic increase in er v with luminosity. 

L05 infer a significant dependence of PVD on color as well 
as luminosity, with red galaxies having a higher dispersions. 
L05 conclude that the dependence of the PVD with luminosity 
results from a population of faint red galaxies in clusters. The 
large difference in the PVD of faint red and blue galaxies can 
also be accounted for through halo occupation models. The 
HOD analysis of SDSS correlation functions in Zeha vi et alJ 
( 2005 ) includes analysis of red and blue galaxies. The derived 
occupation functions predict a higher satellite fraction for red 
galaxies and a larger number of faint red galaxies in rich sys- 
tems (see their Figure 23). These occupation functions would 
produce a substantially higher PVD for red galaxies. The con- 
clusion of L05 is consistent with our interpretation of these 
data through / sat ; faint satellite galaxies are primarily red ob- 
jects, while faint central galaxies are primarily blue galaxies 
in low-mass halos. 

The PVD is physically interesting quantity, both in terms of 
cosmology and galaxy formation, but it cannot be measured 
directly with precision. It can be inferred from the redshift- 
space correlation function by model fitting, but the results are 
model dependent and may be affected by the assumption of 
scale independence of th e PVD. We have in cluded a quali- 
tative comparison to the Hawkins et al. (2003) configuration- 
space estimate of the PVD for a flux-limited sample of 2dF- 
GRS galaxies, but we have focused most of our attention on 
JB and L05's estimate for volume-limited samples, which use 
the dispersion model to infer the PVD form the redshift-space 
power spectrum. We have tested their implementation of this 
model with our fiducial HOD model, demonstrating that the 
effective scale in configurations space probed by a given k 
scale in Fourier space changes with the luminosity of the sam- 
ple, complicating the comparison to predicted trends. At low 
luminosities, Uk at k = 1 Mpc" 1 h probes the regime where 
one-halo pairs dominate the statistics, while for higher lumi- 
nosities this scale probes the the two-halo regime. Therefore, 
while the luminosity dependence of the PVD in configuration 
space is flat at r = 1 /r'Mpc, crj at k = 1 Mpc -1 h follows a 
strong trend with luminosity, comparable to the dependence 
measured by JB and L05. 

For our chosen cosmology, with Q, m = 0.3 and erg = 0.9, the 
normalization of the predicted PVD is too high, both in con- 
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FIG. 8. — Comparison of the fiducial HOD model to the observational results of JB and L05. Open squares and filled circles show the measurements of JB and 
L05, respectively. The dotted line plots the CLF prediction shown in both papers. The large filled squares are the dispersion model results of the fiducial HOD 
model, averaged over the three k values nearest k = 1 Mpc~' h. Error bars are the error in the mean. The solid line shows the HOD results shifted by a factor of 
0.8 to represent a universe with Q m = 0.2. 



figuration space and in Fourier space. This discrepancy can 
be resolved by lowering f2„, to 0.2, thus reducing halo ve- 
locities by (0.2/0.3) 6 and virial dispersions by (0.2/0.3) 05 . 
Our predicted PVD would also be lower for a lower value of 
o"jj. At large scale, cr v scales roughly linearly with erg, while 
in the one-halo regime cr v ~ erg for models constrained to 
match the same w p (r p ). It is difficult to use the PVD mea- 
surements on their own to constrain cosmological parame- 
ters, however, because the normalization of our predictions 
also depends on the velocity bias of both central and satel- 
lite galaxies. A more general model would allow the satellite 
galaxy dispersion to be cr sat = a v adm, and the central galaxy 
dispersion to be er cen = a vc <Jd m - In principle, both a v and a vc 
could depend on halo mass and galaxy luminosity because of 
dynamical effec ts on the orbits of s ubhalos that host satellite 
galaxies. Ivan den Bosch et ail J2005I) argue that a vc ~ 0.25 
for 2dFGRS groups, but in tests we find that significantly 
higher values of a vc are required to alter our predictions. At 
small scales <r v ~ a v , while at large scales er v , ~ a® . Thus 
the normalization of our prediction can be altered by velocity 
bias. For the f2 m = 0.3, erg = 0.9 model used here, a w 0.65 
is required to match the large-scale amplitude of the 2dFGRS 
measurements in Figure [6] a bias that is stronger than that 
predicted by numerical simulations. 

The scale dependence of the P VD predicted by the HOD is 
in general agreement with that of Hawkins et al. ( 2003); both 
the predictions and observations show a smooth transition be- 



tween the one-halo and two-halo regimes as a function of pro- 
jected separation. The <tj curves in both JB and L05 show an 
increase in cr* at k ~ 1 Mpc" 1 h, but the measurements lack 
the strong two-halo to one-halo transition in our HOD models 
(Figure[7J. Changes in a v and erg alter the PVD differently at 
different scales. The difference in scale dependence may im- 
ply a value of a v < 1, or a higher erg than that assumed here, 
although a higher value of erg would exacerbate the normal- 
ization problem in F igure s|6| and 00 

Recent papers bv lGao et al] J2005I) and lHarker et alJ J2005I) 
have shown that halo formation times correlate with the clus- 
tering strength of halos (at fixed mass), implying that for- 
matio n epoch c orrelates with local environ ment. Papers by 
IWechsler et alJ <2005l) and lZhu etaD ||2Q06) demonstrate that 
the number of subhalos correlates with halo formation epoch. 
The implication for the HOD is that (N)m could also corre- 
late with environment, an effect not taken into account in our 
modeling of w p (r p ). If this correlation is strong, the occupa- 
tion functions derived in Figure ^ will t> e systematically off 
their true values. In terms of pairwise velocity statistics, how- 
ever, these errors would not signific antly c hang e the predic- 
tions. First, the effect of iGao et alJ fc005l) and lHarker et alJ 
(2005) is strong for low-mass halos, where (N sat )' M <C 1. Any 
errors accrued would therefore be localized to the central oc- 
cupation function. Because the halo-halo velocity dispersion 
is insensitive to halo mass for masses less than the non-linear 
mass scale, these types of errors would not alter our predic- 
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tions. Second, assuming that the satellite occupation func- 
tions are systematically off by as much as the difference in 
(Nsat)M between adjacent magnitude bins, the change in a v per 
half magnitude is ^ 5% or less for most separatio ns. Lastly, 
tests u sing a hydrodynamic simulation presented in Y oo et alJ 
(2005) demonstrate that any environmental variation of the 
HOD has little impact on the galaxy correlation function and 
is not likely to significantly alter the derived occupation func- 
tions for luminosity-defined samples. The recent numerical 
results do underscore the need to quantify the effect that envi- 
ronmental correlations have on predictions of halo occupation 
models, but we do not expect changes to the qualitative con- 
clusions in this paper. 

Ultimately, there is more information about the PVD in 
the full, two-dimensional redshift-space correlation function 
than in a v from the stre aming model or o> from the disper- 
sion model. Tinker et al. ( 2005a) show that modeling redshift- 
space distortions at small and large scales can separately con- 
strain a v , as, and Q m . Detailed observational comparisons, 
imposing the constraint that erg and fi m must be the same for 



all samples, could reveal luminosity and mass dependence of 
a v , providing physical insight into the dynamical evolution 
of galaxies in groups and clusters. We will apply these ap- 
proaches to 2dFGRS and SDSS observations in future work. 



JT would like to thank Shaun Cole, Charlie Conroy, Carlos 
Frenk, Andrey Kravtsov, Hiranya Peiris, Risa Wechsler, Ren- 
bin Yan, Andrew Zentner, and Zheng Zheng for useful dis- 
cussions. JT would also like to thank the generous hospitality 
of the Institute for Computational Cosmology at the Univer- 
sity of Durham, where part of this work was completed. PN 
acknowledges support of a PPARC PDRA fellowship and the 
computer resources of ETH Zurich. DW acknowledges the 
support of NSF grant AST-0407125. Portions of this work 
were performed under the auspices of the U.S. Dept. of En- 
ergy, and supported by its contract #W-7405-ENG-36 to Los 
Alamos National Laboratory. Computational resources were 
provided by the LANL open supercomputing initiative. 



REFERENCES 



Bahcall, N. A., Dong, E, Bode, P., Kim, R., Annis, J., McKay, T. A., Hansen, 
S., Schroeder, J., Gunn, J., Ostriker, J. R, Postman, M., Nichol, R. C, 
Miller, C, Goto, T., Brinkmann, J., Knapp, G. R., Lamb, D. O., Schneider, 
D. P., Vogeley, M. S., & York, D. G. 2003, ApJ, 585, 182 

Baugh, C. M., Croton, D. J., Gaztanaga, E., Norberg, P., Colless, M., Baldry, 
I. K., Bland-Hawthorn, J., Bridges, T., Cannon, R., Cole, S., Collins, C, 
Couch, W., Dalton, G., De Propris, R., Driver, S. P., Efstathiou, G., Ellis, 
R. S., Frenk, C. S., Glazebrook, K., Jackson, C, Lahav, O., Lewis, I., 
Lumsden, S., Maddox, S., Madgwick, D., Peacock, J. A., Peterson, B. A., 
Sutherland, W., & Taylor, K. 2004, MNRAS, 351, L44 

Bean, A. J., Ellis, R. S., Shanks, T., Efstathiou, G., & Peterson, B. A. 1983, 
MNRAS, 205, 605 

Benson, A. J., Cole, S., Frenk, C. S., Baugh, C. M., & Lacey, C. G. 2000, 

MNRAS, 311,793 
Berlind, A. A. & Weinberg, D. H. 2002, ApJ, 575, 587 
Berlind, A. A., Weinberg, D. H., Benson, A. J., Baugh, C. M., Cole, S., Dave, 

R., Frenk, C. S., Jenkins, A., Katz, N., & Lacey, C. G. 2003, ApJ, 593, 1 
Bullock, J. S., Kolatt, T. S., Sigad, Y., Somerville, R. S., Kravtsov, A. V., 

Klypin, A. A., Primack, J. R., & Dekel, A. 2001, MNRAS, 321, 559 
Cole, S., Fisher, K. B., & Weinberg, D. H. 1994, MNRAS, 267, 785 
— . 1995, MNRAS, 275, 515 

Cole, S., Percival, W. J., Peacock, J. A., Norberg, P., Baugh, C. M., Frenk, 
C. S., Baldry, I., Bland-Hawthorn, J., Bridges, T., Cannon, R., Colless, M., 
Collins, C, Couch, W., Cross, N. J. G, Dalton, G, Eke, V. R., De Propris, 
R., Driver, S. P., Efstathiou, G, Ellis, R. S., Glazebrook, K., Jackson, C, 
Jenkins, A., Lahav, O., Lewis, I., Lumsden, S., Maddox, S., Madgwick, D., 
Peterson, B. A., Sutherland, W., & Taylor, K. 2005, MNRAS, 362, 505 

Colless, M., Dalton, G., Maddox, S., Sutherland, W., Norberg, P., Cole, 
S., Bland-Hawthorn, J., Bridges, T, Cannon, R., Collins, C, Couch, W., 
Cross, N., Deeley, K., De Propris, R., Driver, S. P., Efstathiou, G., Ellis, 
R. S., Frenk, C. S., Glazebrook, K., Jackson, C, Lahav, O., Lewis, I., 
Lumsden, S., Madgwick, D., Peacock, J. A., Peterson, B. A., Price, I., 
Seaborne, M., & Taylor, K. 2001, MNRAS, 328, 1039 

Colless, M., Dalton, G., Maddox, S., Sutherland, W., Norberg, P., Cole, 
S., Bland-Hawthorn, J., Bridges, T, Cannon, R., Collins, C, Couch, W., 
Cross, N., Deeley, K, de Propris, R., Driver, S. P., Efstathiou, G., Ellis, 
R. S., Frenk, C. S., Glazebrook, K., Jackson, C, Lahav, O., Lewis, I., 
Lumsden, S., Madgwick, D., Peacock, J. A., Peterson, B. A., Price, I., 
Seaborne, M., & Taylor, K. 2003, VizieR Online Data Catalog, 7226, 

Cooray, A. & Sheth, R. 2002, Phys. Rep., 372, 1 

Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 
371 

Davis, M. & Peebles, P. J. E. 1983, ApJ, 267, 465 
Doran, M. & Mueller, C. M. 2004, JCAP, 0409, 003 

Dunkley, J., Bucher, M., Ferreira, P. G., Moodley, K., & Skordis, C. 2005, 
MNRAS, 356, 925 

Efstathiou, G., Frenk, C. S., White, S. D. M., & Davis, M. 1988, MNRAS, 
235, 715 

Faltenbacher, A., Kravtsov, A. V., Nagai, D., & Gottlober, S. 2005, MNRAS, 
358, 139 

Gao, L., Springel, V., & White, S. D. M. 2005, MNRAS, 363, L66 
Hamilton, A. J. S. 1993, ApJ, 417, 19 

Harker, G, Cole, S., Helly, J., Frenk, C. S., & Jenkins, A. 2005, MNRAS, 

Submitted (astro-ph/05 10488) 
Hatton, S. & Cole, S. 1998, MNRAS, 296, 10 
— . 1999, MNRAS, 310, 1137 



Hawkins, E., Maddox, S., Cole, S., Lahav, O., Madgwick, D. S., Norberg, P., 
Peacock, J. A., Baldry, I. K, Baugh, C. M., Bland-Hawthorn, J., Bridges, 
T., Cannon, R., Colless, M., Collins, C, Couch, W., Dalton, G., De Propris, 
R., Driver, S. P., Efstathiou, G., Ellis, R. S., Frenk, C. S., Glazebrook, 
K, Jackson, C, Jones, B., Lewis, I., Lumsden, S., Percival, W., Peterson, 
B. A., Sutherland, W., & Taylor, K. 2003, MNRAS, 346, 78 

Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, S., Evrard, 
A. E., Couchman, H. M. P., & Yoshida, N. 2001, MNRAS, 321, 372 

Jing, Y. P. & Borner, G. 2001, MNRAS, 325, 1389 

— . 2004, ApJ, 617, 782 

Jing, Y. P., Mo, H. J., & Boerner, G. 1998, ApJ, 494, 1 
Kaiser, N. 1987, MNRAS, 227, 1 

Kauffmann, G, Nusser, A., & Steinmetz, M. 1997, MNRAS, 286, 795 
Klypin, A. A., Gottlober, S., Kravtsov, A. V., & Khokhlov, A. M. 1999, ApJ, 
516,530 

Kravtsov, A. V., Berlind, A. A., Wechsler, R. H., Klypin, A. A., Gottlober, 

S., Allgood, B., & Primack, J. R. 2004, ApJ, 609, 35 
Landy, S. D. & Szalay, A. S. 1993, ApJ, 412, 64 

Li, C, Jing, Y. P., Kauffmann, G, Boerner, G, White, S. D. M., & Cheng, 

F. Z. 2005, MNRAS, accepted, (astro-ph/0509874) 
Ma, C.-P. & Fry, J. N. 2000, ApJ, 543, 503 

Mandelbaum, R, Seljak, U., Kauffmann, G., Hirata, C. M., & Brinkmann, J. 
2005, MNRAS, submitted, (astro-ph/05 11 164) 

Mohayaee, R. & Tully, R. B. 2005, ApJ, accepted (astro-ph/0509313) 

Navarro, J., Frenk, C, & White, S. 1997, ApJ, 490, 493 

Norberg, P., Baugh, C. M., Hawkins, E., Maddox, S., Madgwick, D., Lahav, 
O., Cole, S., Frenk, C. S., Baldry, I., Bland-Hawthorn, J., Bridges, T., 
Cannon, R., Colless, M., Collins, C, Couch, W., Dalton, G., De Propris, 
R., Driver, S. P., Efstathiou, G., Ellis, R. S., Glazebrook, K., Jackson, C, 
Lewis, I., Lumsden, S., Peacock, J. A., Peterson, B. A., Sutherland, W., & 
Taylor, K. 2002a, MNRAS, 332, 827 

Norberg, P., Baugh, C. M., Hawkins, E., Maddox, S., Peacock, J. A., Cole, 
S., Frenk, C. S., Bland-Hawthorn, J., Bridges, T., Cannon, R., Colless, M., 
Collins, C, Couch, W., Dalton, G., De Propris, R., Driver, S. P., Efstathiou, 

G. , Ellis, R. S., Glazebrook, K., Jackson, C, Lahav, O., Lewis, I., Lumsden, 
S., Madgwick, D., Peterson, B. A., Sutherland, W., & Taylor, K. 2001, 
MNRAS, 328, 64 

Norberg, P., Cole, S., Baugh, C. M., Frenk, C. S., Baldry, I., Bland-Hawthorn, 
J., Bridges, T., Cannon, R., Colless, M., Collins, C, Couch, W., Cross, 
N. J. G., Dalton, G., De Propris, R., Driver, S. P., Efstathiou, G., Ellis, 
R. S., Glazebrook, K, Jackson, C, Lahav, O., Lewis, I., Lumsden, S., 
Maddox, S., Madgwick, D., Peacock, J. A., Peterson, B. A., Sutherland, 
W., & Taylor, K. 2002b, MNRAS, 336, 907 

Peacock, J. A. & Dodds, S. J. 1994, MNRAS, 267, 1020 

Peacock, J. A. & Smith, R. E. 2000, MNRAS, 318, 1144 

Peebles, P. J. E. 1976, ApJ, 205, L109 

Porciani, C. & Giavalisco, M. 2002, ApJ, 565, 24 

Porciani, C. & Norberg, P. 2006, MNRAS, submitted 

Sanchez, A. G. et al. 2005, MNRAS, in press (astro-ph/0507583) 

Sargent, W. L. W. & Turner, E. L. 1977, ApJ, 212, L3 

Scoccimarro, R. 2004, Phys. Rev. D, 70, 083007 

Scoccimarro, R., Sheth, R. K., Hui, L., & Jain, B. 2001, ApJ, 546, 20 

Seljak, U. 2000, MNRAS, 3 1 8, 203 

Seljak, U. & Warren, M. S. 2004, MNRAS, 355, 129 

Seljak, U. & Zaldarriaga, M. 1996, ApJ, 469, 437 

Sheth, R. K. 1996, MNRAS, 279, 1310 

Sheth, R. K., Mo, H. J., & Tormen, G. 2001, MNRAS, 323, 1 

Sheth, R. K. & Tormen, G. 1999, MNRAS, 308, 119 



15 



Slosar, A., Seljak, U., & Tasitsiomi, A. 2006, MNRAS, 366, 1455 

Tinker, J. L., Weinberg, D. H., & Zheng, Z. 2005a, MNRAS, submitted 

(astro-ph/0501029) 
Tinker, J. L., Weinberg, D. H., Zheng, Z., & Zehavi, I. 2005b, ApJ, 631, 41 
van den Bosch, F. C, Mo, H. J., & Yang, X. 2003, MNRAS, 345, 923 
van den Bosch, F. C. et al. 2005, MNRAS, submitted (astro-ph/0502466) 
Warren, M. S., Abazajian, K., Holz, D. E., & Teodoro, L. 2005, ApJL, 

submitted (astro-ph/0506395) 
Warren, M. S. & Salmon, J. K. 1993, in Supercomputing '93 
Wechsler, R. H., Zentner, A. R., Bullock, J. S., & Kravtsov, A. V. 2005, ApJ, 

submitted, (astro-ph/05 1 24 1 6) 
Yang, X., Mo, H. J., & van den Bosch, F. C. 2003, MNRAS, 339, 1057 
Yoo, J., Tinker, J. L., Weinberg, D. H., Zheng, Z., Katz, N., & Dave, R. 2005, 

ApJ, submitted, (astro-ph/05 1 1 5 80) 
Yoshikawa, K., Jing, Y. P., & Borner, G. 2003, ApJ, 590, 654 
Zehavi, I., Blanton, M. R., Frieman, J. A., Weinberg, D. H., Mo, H. J., Strauss, 

M. A., et al. 2002, ApJ, 571, 172 
Zehavi, I., Weinberg, D. H., Zheng, Z., Berlind, A. A., Frieman, J. A., 

Scoccimarro, R., Sheth, R. K., Blanton, M. R., Tegmark, M., Mo, H. J., 

et al. 2004, ApJ, 608, 16 



Zehavi, I., Zheng, Z., Weinberg, D. H., Frieman, J. A., Berlind, A. A., 
Blanton, M. R., Scoccimarro, R., Sheth, R. K., Strauss, M. A., Kayo, 
I., Suto, Y, Fukugita, M., Nakamura, O., Bahcall, N. A., Brinkmann, 
J., Gunn, J. E., Hennessy, G. S., Ivezic, Z., Knapp, G. R., Loveday, J., 
Meiksin, A., Schlegel, D. J., Schneider, D. R, Szapudi, I., Tegmark, M., 
Vogeley, M. S., & York, D. G. 2005, ApJ, 630, 1 

Zheng, Z. 2004, ApJ, 610, 61 

Zheng, Z., Berlind, A. A., Weinberg, D. H., Benson, A. J., Baugh, C. M., 
Cole, S., Dave, R., Frenk, C. S., Katz, N., & Lacey, C. G. 2005, ApJ, 633, 
791 

Zheng, Z., Tinker, J. L., Weinberg, D. H., & Berlind, A. A. 2002, ApJ, 575, 
617 

Zhu, G, Zheng, Z., Lin, W. P., Jing, Y. P., Kang, X., & Gao, L. 2006, ApJ, 

submitted, (astro-ph/0601 120) 
Zurek, W. H., Quinn, P. J., Salmon, J. K, & Warren, M. S. 1994, ApJ, 431, 

559 



